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ABSTRACT 

A  deficiency  is  pointed  out  in  Stevenson's  method  of  reducing  the  solution  of 
electromagnetic  scattering  problems  to  a  succession  of  standard  potential  problems 
whose  solutions  determine  terms  in  the  low  frequency  expansion  of  the  scattered 
field.  An  alternate  approach  is  presented,  for  perfectly  conducting  scatterers, 
which  not  only  removes  the  difficulty  but  also  is  simpler  and  more  explicit  than 
Stevenson's  method.  The  details  of  the  analogous,  though  simpler,  scalar  scatter¬ 
ing  problems  are  also  presented. 


iv 
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I 

INTRODUCTION 

The  purpose  of  this  report  is  to  describe  a  method  of  reducing  scattering 
problems  to  a  series  of  potential  problems .  We  deal  with  a  general  class  of  three 
dimensional  scatterers,  smooth,  closed,  bounded,  in  short  those  surfaces  for  which 
Green's  theorem  in  any  of  its  guises  may  be  invoked.  The  solution  of  a  scattering 
problem,  for  arbitrary  excitation,  is  expressed  as  a  series  in  ascending  powers  of 
wave  number,  k.  This  series  is  known  by  a  variety  of  names,  including  Rayleigh 
series,  quasi-static  series,  and  low  frequency  expansion.  That  the  first  term  in 
such  a  series  could  be  found  as  the  solution  of  a  potential  problem  was  observed  by 
Rayleigh  (1897)  who  determined  this  term  explicitly  for  a  variety  of  scatterers  of 
both  acoustic  and  electromagnetic  waves.  For  scalar  scattering,  the  determination 
of  succeeding  terms  in  this  series  as  solutions  of  potential  problems  has  been  de¬ 
scribed,  in  varying  detail,  by  Noble  (1962),  Morse  and  Feshbach  (1953),  and  Dar¬ 
ling  and  Senior  (1965).  (See  Kleinman  (1965a)  for  a  more  complete  bibliography.) 

The  derivation  of  successive  terms  in  this  series  for  electromagnetic  scat¬ 
tering  was  described  by  Stevenson  (1953a).  Actually  Stevenson  described  two  meth¬ 
ods,  one  for  finding  the  general  term  in  the  series  and  a  second  special  technique 
for  finding  the  first  three  terms.  All  of  his  specific  calculations  (Stevenson,  1953b) 
were  carried  out  using  this  special  technique.  No  attempt  to  utilize  the  general 
method  for  obtaining  higher  order  terms  has,  to  this  writer's  knowledge,  been  re¬ 
ported,  which  indicates  that  if  attempts  were  made,  they  were  unsuccessful.  More 
likely,  there  were  none.  This  is  due  to  the  fact  that  the  analysis  is  sufficiently  in¬ 
volved  to  discoura^s  most  efforts  to  derive  more  than  three  terms  in  a  low  fre¬ 
quency  expansion  (that  Stevenson  treats  the  more  general  case  of  penetrable  scat¬ 
terers  certainly  doesn't  help).  For  these,  Stevenson's  special  simpler  technique 
suffices.  An  attempt  to  clarify  the  Stevenson  method  was  made  by  Senior  and 
Sleator  (1964)  and  the  present  report  may  be  considered  an  outgrowth  of  their  work. 
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The  present  work  demonstrates  that  the  method  proposed  by  Stevenson  for 
finding  the  general  term  in  the  series  needs  clarification  at  best  and  at  worst  leads 
to  incorrect  results.  An  alternate  method,  preset ving  the  spirit  of  Stevenson's 
approach  and  indeed  largely  based  on  it,  is  presented  which  hopefully  embodies  both 
clarity  and  correctness.  Conciseness  has  been  sacrificed  in  an  attempt  to  minimize 
the  chances  of  further  obscuring  the  subject. 

The  procedure  in  the  electromagnetic  (vector)  case  is  a  natural  extension  of 
the  technique  employed  in  the  scalar  case.  For  this  reason,  and  also  to  introduce 
some  notation  as  well  as  concepts  in  the  simplest  setting,  the  next  section  is  de¬ 
voted  to  a  discussion  of  how  scalar  scattering  problems  may  be  reduced  to  the  study 
of  a  succession  of  potential  problems.  In  Section  3  we  describe  Stevenson's  method 
for  treating  the  analogous  vector  problem  and  show  why  it  is  unsatisfactory.  Sec¬ 
tion  4  presents  an  alternative  to  Stevenson's  method  which  eliminates  its  short¬ 
comings.  Section  5  is  devoted  to  an  illustrative  example. 
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II 

SCALAR  SCATTERING 

In  this  section  we  show  how  a  scalar  scattering  problem  with  Dirichlet  or 
Neumann  boundary  conditions  may  be  reduced  to  a  succession  of  "standard"  poten¬ 
tial  problems.  These  terms  will  be  precisely  defined  as  they  are  introduced. 

Let  B  denote  the  boundary  of  a  smooth,  closed,  bounded  surface  in  Euclid¬ 
ian  3-space  (or  the  union  of  a  finite  number  of  such  surfaces  provided  they  are  dis¬ 
joint),  let  n  denote  the  outward  drawn  unit  normal  at  any  point  of  B  and  let  V  be 
the  volume  exterior  to  B.  Erect  a  cartesian  coordinate  system  with  origin  in  B  and 
let  r  denote  a  radius  vector  to  a  general  point  (x,  y,  z)  and  r  denote  a  point  on  B . 

Furthermore  denote  by  R  the  distance  between  r  and  r  ,  i.e., 

B 

R  =  =  " *b)2  +  (y " yB^  +  (z  ‘  ZB)2  ‘  *21) 

The  geometry  is  illustrated  in  Fig.  1. 
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By  a  scalar  scattering  problem  for  the  surface  B  is  meant  the  determination 
of  how  the  presence  of  the  surface  perturbs  an  incident  field,  0^°,  that  is,  finding 
a  function  $(  r )  such  that 


(V2+k2;^  =  0 


reV  , 


r-*a> 


(2.2) 


(2.3) 


and  either 


$<?„>  =  -  $mc<?  j  . 


(2.4a) 


9$<r)  _  3^iaC(r) 


(2.4b) 


Equation  (2.3)  is  a  statement  of  Sommerfeld's  radiation  condition  which  im-  ! 

eikr 

plies  that  outgoing  waves  look  like  — ~  f (0, 0)  for  large  r.  The  boundary  condi¬ 
tions  (2.4a)  and  (2.4b)  are  Dirichlet  and  Neumann  conditions  respectively.  Specify¬ 
ing  either  one  is  sufficient  to  guarantee  the  existence  of  a  unique  function  $  hence 
both  the  values  of  the  function  and  its  normal  derivative  may  not  be  assigned  arbi¬ 
trarily.  We  will  consider  the  Dirichlet  and  Neumann  problems  separately  but  the 
analysis  is  quite  similar. 

The  starting  point  is  the  Helmholtz  integral  representation  of  regular  solu¬ 
tions  of  (2.2);  viz, 


re>  ■  *  Si 


an  R 


ikR  a(K 
R~  an 


(2.5) 


The  integration  is  carried  out  over  the  surface  and  the  normal  derivative  is 
n  •  V.  Next  we  assume  that  the  unknown  function  $  may  be  expressed  as  a  conver- 
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gent  power  series  in  k.  Actually  this  need  not  be  assumed,  that  is,  it  may  be 
proven  that  there  does  exist  such  an  expansion,  convergent  for  k  sufficiently  small 
(see  Werner,  1962  and  Kleinman,  1965b).  It  should  be  noted  that  we  are  consider¬ 
ing  k  real  and  positive  though  the  results  may  be  extended  to  include  complex  val  - 
ues  of  k.  We  write  the  expansion 


(2.6) 


where  the  factor  i  is  included  in  the  expansion  parameter  merely  as  a  convenience. 
The  functions  are  independent  of  k  and  each  of  them  may  be  determined  as 


follows. 


Since  e  is  an  entire  function,  the  series 


oo  . 
ikR  =  ^  (ikR) 


(2.7) 


converges  for  all  k.  Substituting  (2.6)  and  (2.7)  in  (2.  5),  we  obtain 


Jc  Lm=0  t  =0 


t=0  m=0  J  (2. 8) 

As  long  as  Jk  J  is  strictly  less  than  the  radius  of  convergence  we  may  interchange 
summation  and  integration  and  reorder  terms  in  the  double  series  obtaining 

m=0  1=0  m=0  .L  u  — 1  • 


(2.9) 
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Equating  coefficients  of  like  powers  of  k,  with  an  obvious  change  in  notation. 


yields 


*/*  >  -  £  i  irsi  i  R,'m'1  -  i  UdB  •  <2- 10> 

i  =  0, 1,2,  . .. 

In  order  to  determine  (pj  we  must  employ  the  boundary  conditions,  hence  we 
must  distinguish  between  the  two  problems  under  consideration.  Whether  the  inci¬ 
dent  field  is  a  plane  wave,  point  source,  or  linear  combination  of  such  sources  it 
remains  true  that  the  representation  of  the  incident  field  is  analytic  in  k.  Thus  we 
may  write 


uu 

$toc(r>  =  . 


(2.11) 


The  boundary  conditions  (2.4a)  and  (2.4b)  then  imply  that  either 


(2.12a) 


9n  a  .x 
r=r. 


-  -2-  6mC(^) 

8n  '  *  a 

r=r. 


(2.12b) 


Consider  first  the  Dirichlet  problem,  (2, 12a).  Inserting  the  boundary  values  in  the 
integral  representation  (2. 10)  produces  the  system  of  equations 


l  4n  (f-m)'. 


i4ino  a 

8n 


R(-m-l  +  Rl-m 


-1f  $  \ 

9n  -  mj 


(2.13) 
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We  treat  first  the  case  when  i  —  0.  Equation  (2. 13)  becomes  simply 


i(?>  ■  -  \ 

'B 


ac  3  1  J_ 

3n  R  4tt 


(2.14) 


The  unknown  term  on  the  right  is  clearly  an  exterior  potential  function  or  in  the 
language  of  potential  theory  (e.g.  Kellog,  1929)  the  field  of  a  single  layer  distribu¬ 
tion  of  density  3^/dn.  That  is,  if  we  designate  by  0  the  unknown  function. 


J. 

R 


(2.15) 


then  0^  satisfies  the  equation 

V20o<?)  =  0  ,  reV  , 

and  is  regular  at  infinity  in  the  sense  of  Kellog,  viz . 


lim 
r  — >oo 


<  O'  , 


lim 
r— >  oo 


<  oo  . 


(2.16) 


(2.17) 


Furthermore,  with  the  boundary  condition  (2. 12a)  and  the  expression  (2. 14), 


the  values  of  0q  on  B  are  specified,  i.e.. 


W  -  -«X>+ 


B 


inc  _S_ 
3n 


(i) 


dB 


(2.18) 


Note  that  the  integration  in  (2. 18)  must  be  carried  out  before  the  limit  is  taken  so 
the  integrand  is  always  defined.  With  this  proviso  the  right  hand  side  of  (2. 18)  is 
well  behaved  and  completely  specified  in  terms  of  the  incident  field.  Equations 
(2. 16)- (2. 18)  constitute  a  standard  exterior  Dirichlet  potential  problem  which  has  a 
unique  solution.  Next  we  show  that  succeeding  terms  may  be  written  in  terms  of 
solutions  of  similar  problems.  To  this  end  assume  that  (j>Q,  <^,  . .  ,<^  are  all 

»  n  _ _ _ 
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known.  Then  (2. 13)  may  be  written 

fy?>  -  F^rl  +  tyr) 


(2.19) 


where 


Fp) 


_  JL  \  ^  1  \  xinc  8  f-m-1 

4 it  \  ^m  8n 

m-0  JB 


.  V7  1  \  aJL 

4ir  JB  90 


dB 


and 


ft(i>  =  - 


A±ba 

4j  J  H  3n 
JB 


dB 


With  the  assumption  that  all  <$'s  are  know  up  to,  but  not  including  <^,  Ff(r )  is  a 
known  function.  Clearly  0  ( r )  is  again  a  single  layer  distribution,  satisfies  (2. 16) 
and  (2. 17),  and  is  uniquely  determined  with  the  boundary  condition 


tyrB)  =  -$nC(fB,-Fl(?B)  <2'20) 

Again,  care  must  be  used  in  letting  r-*rn  in  one  term  of  F,  but  there  is  no  in- 

b  1 

trinsic  difficulty.  Thus  <£  is  determined  in  terms  of  a  known  function  F  and  a 
solution  of  a  standard  exterior  Dirichlet  potential  problem,  0  .  We  have  shown  that 
is  true  for  i  =0,  and  also  for  i  >0  provided  ^  . .  .<$  have  previously 

been  found.  The  solution  of  the  Dirichlet  scattering  problem  is  then  given  by  (2.6). 

An  exactly  analogous  procedure  may  be  followed  for  the  Neumann  boundary 
condition,  (2. 12b).  Corresponding  to  (2.13)  we  have 


L 
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4* 


1  ^  '  1  \  j"Li-m- 1  9  zinc  .  x  9  _f-m“l"\ _ 

■to  2-gWSi  \  1R  J  “ 

m=0  JB 


Furthermore  this  may  be  written 

$^(r)  =  Gj(r)  +  0^(r) 

where 

^  l 

„i-m-l  6  line 
R  *}„.  ® 


(2 . 21) 


(2.22) 


Z. 

~  1  V  1  \  «i-m-l  9  Xinc 

G<(r)  -  to  A-  (T^mj; 


m  r 

1  \  1  \  r  9  _f-m-l 

4*  (I-  m)'.  \  ^m  9n  R  63  ’ 


the  second  sum  is  identically  zero  if  i  =  0,  and 

*i(?)  =  ir^<^(i)dB  • 


is  completely  determined  U  $1-1'  •••$o  are  known  and  r )  is  a  double  layer 
distribution.  That  is,  ip  is  the  solution  of  a  standard  exterior  Neumann  potential 
problem,  namely. 


V2  djCr)  = 


^(r)  =  0  , 


reV 


ip  regular  in  the  sense  of  Keliog,  (2.17), 


and 


mm 

r=r 


B 
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We  have  thus  demonstrated  that  for  either  Dirichlet  or  Neumann  scattering 
problems  successive  terms  in  the  low  frequency  expansion  may  be  determined  by 
solving  a  succession  of  standard  potential  problems.  That  is,  the  first  term  is  the 
solution  of  such  a  potential  problem  the  second  term  is  expressed  in  terms  of  the 
first  and  the  solution  of  a  potential  problem,  the  third  is  given  in  terms  of  the  first 
two  and  a  potential  solution,  etc . 

Before  closing  this  section,  a  word  should  be  said  about  low  frequency  ex¬ 
pansions  of  the  far  field.  The  Rayleigh  series,  (2.6),  may  be  considered  as  an 
expansion  of  the  near  field  which,  if  all  terms  are  included,  is  also  valid  in  the  far 
field.  If  only  a  finite  number  of  terms  are  known,  tnen  the  truncated  series  does 
not  in  itself  give  much  useful  information  about  the  far  field.  Such  information  is 
available  if  we  again  make  use  of  the  integral  representation  (2.5).  To  this  end 
note  that  for  large  r 


r  •  r 


..  lkr  -  lk - 

lkR  r 

e  e 


B 


ikr  -  ik  r  •  r 


B 


R 


r  =  |r  |,  r  =  r/r  (2.23a) 


and 


ikR 


ikr  -  ik  r  ■  ? 


R 


..  a  e 
'v  xkr  — 


B 


(2.23b) 


If  (r,  Q,  0)  and  (r  ,  6  ,  0  )  are  spherical  coordinates  of  points  r  and  r  respec- 
b  b  b  B 


tively  then 


A  — 

r  •  r 


B 


B  l_ 


fcos8cos0 


B 


sin  9 sin  0  cos(0  -  0  ) 
B  B  J 


(2.24) 


Substituting  (2. '23a)  and  (2.23b)  in  12. 5)  we  obtain,  for  large  r, 


ikr 

e 

47rr 


(2.25) 
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-ikr  •  r. 


—  B  « 

Now  if  we  substitute  expansions  of  <j>  and  e  in  the  right  hand  side  of  (2 . 25)  | 

and  rearrange  terms  we  obtain 

B  (2.26) 

where  $  =  0. 

! 

Examination  of  equation  (2. 26)  reveals  that  knowledge  of  a  finite  number  of  terms  in 

the  low  frequency  expansion  of  the  near  field  (the  d>  's)  provides  similar  informa- 

m 

tion  about  the  low  frequency  expansion  of  the  far  field  coefficient,  i.e.  the  coeffi- 
cient  of  e  /r.  More  specifically,  in  the  Dirichlet  case,  when  the^oundary  con¬ 
ditions  specify  d>  on  B  for  all  m,  then  knowledge  of  the  first  £  ~~~ ’s 
m  9n 

C'%  8<*i  . . 


\dn~  ’  ~dn~  ’  “  *  ~dn  /  wi^  Provide,  with  equation  (2.26),  the  first  t  terms  of 

a§m 

the  far  field  expansion.  In  the  Neumann  case,  (  — —  on  B  given  for  all  m)  the 

an 

first  i  ©  's  will  apparently  give  f+1  terms  in  the  far  field.  However,  it  may  be 
m 

shown  that,  whether  §  is  a  plane  wave  or  a  point  source. 


(2.27) 


hence  only  f  terms  in  the  far  field  are  specified.  Or  another  wray  of  saying  this  is 
that,  in  the  Neumann  problem,  £  near  field  terms  produce  f  + 1  far  field  terms  but 
the  first  term,  i.e.,  the  coefficient  of  (ik)^,  is  always  zero. 
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m 

STEVENSON’S  METHOD 

In  this  section  we  shall  describe  Stevenson’s  attempt  to  generalize  the 
approach  of  Section  2  to  electromagnetic  scattering  and  pay  particular  attention  to 
the  shortcomings,  rather  than  the  strong  points  (which  are  numerous)  of  Stevenson’s 
work.  To  effect  some  simplification,  we  shall  treat  only  the  case  of  scattering  by  a 
perfectly  conducting  surface  whereas  Stevenson  considered  more  general  scatterers. 

1 

It  seems  clear,  however,  that  both  the  criticism  in  this  section  and  the  correction 
in  the  following  section  may  be  applied  in  the  more  general  case. 

The  surface  geometry  and  notation  are  the  same  as  introduced  in  Section  2 
and  depicted  in  Fig.  1,  which  is  here  reproduced  for  convenience. 


By  an  electromagnetic  scattering  problem  for  the  perfectly  conducting  sur¬ 
face  B  is  meant  the  problem  of  determining  how  the  presence  of  the  surface  per- 

inc  Hnc 

turbs  an  incident  electromagnetic  field,  (E  ,  H  ).  That  is,  we  seek  a  solution 
of  Maxwell's  equations 

VxE  =  ikH,  V-  E  =  0  ,  re  V  (3.1) 

VxH  =  -ikE  ,  V  •  H  -  0 
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subject  to  the  boundary  conditions 


nxE 


L  .  =-nxEtac 

rr=rB 


|  Tar 


A  ,, 

n  •  H 


B 


r=r 


a  ^inc 
=  -n  •  H 


B 


r=r. 


(3.2) 


B 


and  the  radiation  condition 
lim  j. 

rx(VxE)  lkr  E 


ao 


lim  a  .  — ,  ^ 

rx(VxH)+ikrH  =  0 

r— *ao 

uniformly  in  r 


(3.3) 


(The  divergence  conditions  and  boundary  condition  on  H  are  redundant.,  i.e.  may 
be  deduced  from  the  other  conditions . ) 

In  attempting  to  show  how  to  reduce  this  problem  to  that  of  solving  a  series 
of  potential  problems,  the  procedure  parallels  that  followed  in  the  scalar  problem. 
Corresponding  to  the  Helmholtz  integral  representation  (2.2)  we  employ  the  expres¬ 
sion  derived  by  Stratton  and  Chu  (see  Stratton,  1941)  which  expresses  the  field  at 
any  exterior  point  in  terms  of  its  values  on  the  surface  B.  [wilcox  (1956)  also  de¬ 
rives  these  formulas  but  strangely  omits  any  reference  to  the  Stratton-Chu  workj 


E(r)  * 


ikR  .. 

e _  a  •  *  ,  lk 

— ~  nxEdB  +  — 
R  4?r 


ikl* 

e  a  *  1  _ 

— — ~  nxHdB-— V 
R  47r 


ikR 

> 

~R 


n-  EdB 


(3.4a) 


H(r)  = 


_  nxHdB 
R  47r 


A  '  TS 


4  R 


nxEdB  -  ~  V 
47T 


ikR 

e  a  -* 

—  n-HdB. 

(3.4b) 


Recall  that  R  is  a  function  of  the  coordinates  of  two  points  r  and  ?  ,  everything 
else  in  the  integrands  on  the  right  hand  sides  is  a  function  of  the  integration  vari¬ 
ables  (coordinates  of  r  )  and  V  operate*:  on  r.  For  future  use,  we  denote  bv  V 

B  J  B 

the  operator  onr.  and  note  that 
15 
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ikR  ikR 

V- -  =  _  y  - - 

R  B  R 

Now  following  Stevenson  as  well  as  the  procedure  in  the  scalar  case  we  assume  that 
•«*» 

E  and  H  may  be  expanded  in  series  of  powers  of  k,  i.e.. 


CD  oo 

E(r)  =  y^E  ( r)(ik)m,  H(r)  =  >  »  (r)(ik)m  . 

na  m 


(3.6) 


m=0 


m=0 


As  before,  tois  assumption  has  been  proven  (Werner,  1963),  that  is,  it  is  no  longer 
an  assumption  but  a  consequence  of  (3.1),  (3.2)  and  (3.3).  It  is  perhaps  worthy  of 
note  that  the  reason  this  entire  discussion  concerns  three-dimensional  scattering 
problems  is  that  convergent  expansions  of  the  form  (3. 6)  do  not  exist  for  two-dimen 
sional  scattered  fields . 

Next  we  expand  the  free  space  Green's  function,  e  /R,  in  a  series,  vi z, 


ikR 

ft 

~R~ 


oo 

Y  (ik)  R 
Z — i 
t^o 


(3.7) 


then  substitute  (3.6)  and  (3.7)  in  (3.4a),  (3.4b).  After  interchanging  summation  andj 
integration,  reordering  terms  and  equating  like  powers  of  (ik),  we  obtain 


E/<?> 


4 v  *■ — t1  mi 
m=0 


nxE  Rm_1dB 
i-m 


*  t,  T,ttt-1  __ 

nxH.  R  dB  . 

i-m-1  4tt  ‘ — rf  mi 

m=0 


-  “  v\  n  •  E  Rm  l 

t7r  c — e  ml  1  i-m 


dB! 


I  =  0,1,2,  .. 


(3.8a) 
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H  (r)  =  \  nxH  Rm_1dB-  ~  \  nxE  Rm_1dB 

t  4tt  — r*  m'.  i-ui  47r  r*  m'.  \  i-m-1 

m=0  JB  m=0  JB 


JL 

47 r 


y> 


m=0 


V  \  n-  H  Rm_1  dB 
ml  1  i-m 


;B 


(3.8b) 


l-l 


where  /  -  0  when  1  =  0. 

i=0 

Furthermore,  substituting  the  series  (3.6)  in  Maxwell’s  equations  (3. 1)  yields 


VxE  =  0 
o 

V  •  E  =  0 


VxH  =  0 
o 

VxiIi=-SM 

V  •  H  =  0 
1 


1  =1,2,3... 

1  =  0,1,  2, 3 


(3.9) 

(3.10) 

(3.11) 


and  the  boundary  conditions  (3.2)  become 


nxE, 


=  -fixE^I 

n-  H 

\  *  uincl 

=  -n  •  H 

-  ^  1 
r~rB 

A 

r'rB 

1 

15-t 

ll 

•1»H 

(3.12) 


r=r 


B 


These  last  equations  result  from  the  fact  that,  as  with  scalar  sources,  represen¬ 
tations  of  electromagnetic  plane  waves  or  point  sources  are  analytic  in  k.  There  is 
one  more  condition  of  importance.  With  Maxwell's  equations  and  Stokes' theorem  it 
is  a  simple  matter  to  show  that 

r 


n  •  E  dB  =  0  , 


n-  HdB  =  0 


(3.13) 


'B  °B 

It  follows  then  from  the  series  expansions  (3.6)  that 


n  •  EidB  =  0  , 


n  •  dB  =  0 


(3.14) 


IB 


JB 


f  =  0,  1,2,3, 
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Stevenson  then  proceeds  to  show  how  the  zeroth  order  terms,  E  and  H  , 

o  o 

may  be  determined  as  solutions  of  potential  problems .  For  the  perfectly  conducting 

case  this  reduction  to  potential  problems  for  the  zeroth  order  terms  will  be  included 

in  the  general  treatment  of  the  following  section,  and,  since  we  have  no  quarrel  with 

Stevenson's  results  for  these  terms,  the  details  will  be  omitted  here.  To  calculate 

higher  order  terms,  Stevenson  proposes  the  following  procedure:  Suppose 

Eq,  Ej,  H^,  are  known.  To  find  or  H^,  determine  first 

a  particular  solution  of  (3.10),  that  is,  find  functions  F  and  G  such  that 

f  t 

ViV“m  3,1,3  Vx2i  =  -e<.  1 


The  differences  between  these  particular  solutions  and  the  true  coefficients, 
E^-  and  -  G^,  are  gradients  of  unknown  potential  functions  (not  necessarily 
regular  at  infinity)  i.e., 

Ej  =  (3.16a) 

=  (3.16b) 

Substitute  (3.16a)  and  (3.16b)  into  the  integral  expressions  (3.8a)  and  (3.8b)  respec¬ 
tively,  also  introduce  the  boundary  conditions,  (3.12).  There  results  equations  for 
E^  and  which  contain  some  known  terms  and  some  unknown.  It  is  then  possible 
to  show  that  the  unknown  terms  are  now  exterior  potential  functions  (regular  at  in¬ 
finity)  which  may  be  determined  as  solutions  of  standard  potential  problems. 

The  process,  once  begun,  appears  to  be  both  correct  and,  in  the  details  of 
its  execution,  ingenious.  The  source  of  trouble,  however,  is  right  at  the  beginning; 
namely  how  does  one  determine  particular  solutions  of  the  equations 


VxE 


l 


VxH 

1 
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In  a  separate  paper,  Stevenson  (1954)  points  out  that  necessary  and  sufficient! 
conditions  for  the  equation 


VxF(r)  =  f(r)  , 


reV 


to  have  a  solution  are 


V-  f  =  0 


and 


n  '  f  ds  =  0 


(3.17) 


(3.18) 


(3.19) 


'B 


where,  if  B  consists  of  a  number  of  disjoint  surfaces,  B^, . .  .B.,  then  (3.19)  must 
hold  for  each  separately,  as  well  as  the  sum.  With  this  we  have  no  quarrel.  Stev¬ 
enson  then  goes  on  to  assert  that  an  explicit  solution  of  the  problem  is  given  by 


F(r)  =7-  Vx 
47 r 


'all 

space 


f(r  ) 

v 

.  dv 
K(r,  r  ) 
v 


(3.20) 


provided  that  f  satisfies  (3.18)  and  (3.19). 

Since  the  integration  is  over  all  space,  not  merely  V,  the  exterior  of  B, 
this  expression  requires  some  explanation.  In  the  first  place,  f(r)  is  originally 
defined  only  exterior  to  B.  To  extend  the  definition  to  the  interior,  Stevenson  pro- 
poses  to  choose  f  so  that  (3. 18)  remains  true  and  that  n  •  f  is  continuous  at  B. 

This  he  accomplishes  by  choosing 


f(r)  =  Vu(r)  , 


where 


V2  u  =  0  , 


r  interior  to  B 


r  in  B 


(3.21) 


(3.22) 
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n  •  Vu 


r=r 


3u 

3n 


B 


n.  f 


r=r. 


r=r 


B 


(3.23) 


This  is  a  standard  inferior  Neumann  potential  problem  for  u  and  has  a  unique  solu¬ 
tion  provided  that 


(3.24) 


That  (3.24)  holds  is  guaranteed  by  (3. 19).  Thus  the  extension  to  the  interior  is 
carried  out,  once  tnis  potential  problem  is  solved.  Equation  (3.20)  then  is  the  re¬ 
quired  solution  provided  the  integral  exists,  that  is,  provided 

_k  3 

f  =  0(1/ r  )  as  r->oo  .  (3.25) 

Stevenson  describes  the  proof  and  we  shall  demonstrate  it  in  detail  in  the  following 
section  where  we  again  make  use  of  this  device.  Now  however,  we  accept  it  and 
finally  get  to  the  heart  of  the  matter,  namely,  what  do  we  do  if  f  is  defined  origin¬ 
ally  in  the  infinite  region  V,  but  does  not  satisfy  the  necessary  order  condition  at 

infinity,  equation  (3.25)?  This  in  fact  is  exactly  what  happens  since  E  and  it 
3  if 

vanish  as  l/r  only  for  i  =  0  which  allows  us,  using  the  method  described,  to  de- 

—A 

termine  E  and  H  but  apparently  no  higher  order  terms.  (Actually  we  may  go  one 
1  *  2 

term  further  since  the  l/r  terms  don't  contribute  to  the  integral.)  Stevenson  was 
aware  of  this  and  proposed  the  following  procedure: 

If  V  is  the  unbounded  region  exterior  to  B  and  if  f  does  not  vanish  at  infin¬ 
ity  to  the  required  order,  first  surround  B  by  a  surface  B  .  Then  redefine  f  ex- 

o 

terior  to  B^  in  terms  of  the  solution  of  an  exterior  potential  problem,  namely,  let 

f(r)  =  Vu,  r  exterior  to  B  (3.26) 

o 
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where 


V2  u  —  0 


3u 

On 


/*  r 
n  •  f 


r=r 


B, 


•i  -jk 

r=r 


(3.27) 


B 


o 


u  regular  at  infinity. 

This  problem  has  a  unique  solution  u  and,  since 


r\ 

n-fdB  =  \  ~dB  =  0 


dn 


(3.28) 


'B 


'B 


it  follows  that  u  =  0(i/r  )  hence  f  will  satisfy  (3.25).  With  f  thus  redefined,  the 

solution  (3.20)  exists  and  is  valid  in  the  portion  of  V  interior  to  B  where  B  can 

o  o 

be  taken  arbitrarily  large. 

With  that,  Stevenson  apparently  considers  the  subject  closed.  The  implica¬ 
tion  is  that  since  B^  may  be  taken  arbitrarily  large  we  may  take  it  as  a  sphere 
whose  radius  becomes  infinite  and  then  (3.20)  will  represent  the  solution  we  seek 
throughout  V.  But,  unfortunately,  if  f  were  a  function  whose  original  behavior  at 
infinity  was  insufficient  to  guarantee  existence  of  the  integral  in  (3.20),  then  the 

limit  of  the  integral  with  f  redefined  may  not  exist  as  the  radius  of  B  becomes  in- 

- -  Q 

finite.  This  argument  by  which  the  unpleasant  behavior  at  infinity  is  avoided  (that 
is,  confining  attention  to  a  finite  volume,  carrying  out  the  calculation,  and  then 
letting  the  volume  become  infinite)  is  not  only  employed  by  Stevenson  but  others  as 
well,  e.g.  Morse  and  Feshbach  (1953,  I,  p.  53).  It  does  produce  the  desired  results 
in  many  cases.  For  example,  the  process  is  valid  whenever  f  is  the  gradient  of  a 
potential  function,  regardless  of  its  behavior  at  infinity  (which  includes  the  example 


:,tIf  B^  is  any  surface  entirely  containing  B  and  equations  (3. 18)  and  (3.19)  hold, 
then  (3.28)  follows  from  Gauss'  theorem  relating  volume  and  surface  integrals. 
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used  by  Morse  and  Feshbach).  That  it  may  also  yield  unacceptable  results  is  illus¬ 
trated  in  the  following  example,  where  f  is  not  ibe  gradient  of  a  potential  function 
though  still  satisfies  (3.18)  and  (3.19).  This  indeed  is  representative  of  the  be¬ 
havior  one  would  encounter  in  actually  attempting  to  find  particular  solutions  of 
(3.10). 


_  m*  m-2 ,  a  *  . 

i\r)  =  Vxr  l  =  mr  (zi  -yi  )  , 
x  y  z 


(3.29) 


f  2  ,  2  2 

[X  +y  +z 


Clearly  a  particular  solution  of  V  x  F  =  f  is 

F  =  rmi 


(3.30) 


However,  let  us  attempt  to  determine  a  particular  solution  using  equation  (3.20). 
First  of  all  it  is  a  trivial  calculation  to  observe  that  (3. 18)  and  (3.19)  are  satisfied 
with  this  particular  f.  For  this  simple  example  we  have  no  scattering  surface  B, 
but  with  Gauss’  theorem  it  is  clear  that  for  any  closed  surface  B,  ^n  •  fdB  =  0. 

Furthermore,  the  function  f  clearly  misbehaves  at  infinity  so  that  to  use  (3 . 20)  we 

must  employ  the  redefinition  of  f.  Thus  choose  B^  to  be  a  large  sphere  of  radius 

r  .  Next  define 
o 


f  =  f 


r  <r 


=  Vu 


r  ^  r 


where 


V  u  =  0 


n  •  Vu 


r  >  r 


*  7 

=  n  •  f 


u  regular  at  infinity. 
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n  „  m-1 

^ 2  m*  .  m  +  2  m*  mxr  a 

F(r)  =  -  -  r  i  +  — —  r  1  - - —r—  r 

3  o  x  m+3  x  m+  3 


r  <  r 


(3.37) 


3(m  +  3)r 


t  <i  -3  —  r) 
3  x  r 


r  >  r 


It  is  a  simple  calculation  to  show  that  this  F  is  indeed  a  solution  of  VxF  =  f  when 

r  <  r  .  However,  it  is  also  clear  that  F  as  defined  in  (3. 37)  does  not  exist  as 
o 

r  — ►cd.  Furthermore,  if  r  remains  finite,  then  the  function  F  not  only  exhibits 
o  o 

an  unwanted  dependence  on  an  arbitrary  parameter  (  the  radius  r  of  B  )  but  also 

o  o 

is  discontinuous  on  B^.  This  violates  the  tacit  requirement  that  F  be  a  differen¬ 
tiable  solution  of  VxF  for  aU  points  in  V. 

How  then  do  we  proceed  in  those  cases  when  Stevenson's  scheme  for  finding 
particular  solutions  apparently  fails?  One  method  would  be  to  attempt  to  show  that 
the  undesirable  part  of  F  is  the  gradient  of  a  scalar  function  and  can  therefore  be 
neglected;  the  remaining  part  of  F  would  still  be  a  solution  of  VxF  =  f.  In  the 
example  above  it  is  easily  seen  that  F  may  be  written 


1.0  m-l 

*  m  +  2  Ida  mxr  a  _ 

F  =  — —  r  l  -  - — —  r  +  V 

m+3  x  m+3 


(-•>?■) 


(3.38) 


hence  a  particular  solution  of  Vx  F  =  f,  valid  everywhere  in  V  may  be  obtained 
merely  by  deleting  the  term  V  |  r™x  j  .  In  general,  however,  the  process  of 
identifying  the  unwanted  terms  with  the  gradient  of  some  functions  may  not  be  so 
easily  accomplished  and  in  any  event  adds  yet  another  complication  to  an  already 
involved  procedure . 

Rather  than  attempt  to  prove  that  this  procedure  can  be  made  correct  in  the 
manner  indicated,  we  shall  end  this  section  having  demonstrated  that,  as  it  stands, 
Stevenson's  procedure  is  ambiguous.  In  the  next  section  we  shall  show  that  this 
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problem  of  finding  particular  solutions  of  (3.10)  may  be  avoided  entirely  and  the 
process  of  determining  successive  terms  in  the  low  frequency  expansion  may  be 
made  more  straightforward. 
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rv 

ELECTROMAGNETIC  SCATTERING  -  AN  ALTERNATE  APPROACH 

In  this  section  we  again  treat  the  problem  of  extending  to  electromagnetic 
scattering  the  method  of  Section  2  whereby  scalar  scattering  problems  are  reduced 
to  a  series  of  standard  potential  problems.  Though  the  method  described  here  de¬ 
parts  from  Stevenson's  approach,  the  debt  to  his  work,  both  in  ideas  and  technique, 
is  large. 

We  formulate  the  problem  exactly  as  in  Section  3  and  the  details  will  not  be 
repeated.  The  starting  point  for  this  analysis  is  the  integral  representation  of  the 
coefficients  in  the  low  frequency  expansions  of  the  scattered  field,  equations  (3.8a) 
and  (3.8b).  That  is,  we  write  the  field  scattered  from  the  perfectly  conducting  body 
B,  as 


oo 

E(r)  =  y^ElrHik/  , 
1=0 


co 

H(  r )  =  /  .iUrW, 

i=0  1 


(4.1) 


then  the  boundary  conditions  at  the  surface  B  are 


*  £  £  £inc 

nxE,  =  -nxE,  , 

l  l  ' 


A  it  A  ulnC 
n  •  H„  =  -n  •  H„ 
l  i 


l  =  0, 1,2. 


(4.2) 


and,  furthermore. 


n  •  EfdB  =  0 
'B 


n-  H.dB  =  0 


t  =  0, 1,2,  . ..  .  (4.3) 


B 


The  Stratton-Chu  integral  representation,  after  expanding  in  powers  of  k, 
equating  coefficients  and  using  the  boundary  conditions  may  be  written 


E,(r)  =  F^(r)  - 


R 


dB 


(4.4a) 


L 


24 


THE  UNIVERSITY  OF  MICHIGAN 

7133-4-T 


nxH, 


H^(  r )  =  G^(  r*)  +  ~  V  x  \  -j^dB 


(4.4b) 


'B 


where 


Fix)  =  -  ~ 
I  4tt 


“*®  m-  JB 


inc  „m-l 
R  dB 
m 


(4. 5a) 


i-i  c  if 

+  7-  /  .  A  \  nxH  Rm'1dB-  -r-^^7  v\  n-E  Rm'' 
m=0  JB  m=l  JB 


dB 


and 


1 


l 

c - 1  m ' 


m=l 


(4.5b) 


1-1 


-y- 

4 TT  — r1  m'. 
m=0 


JL 

a  =*inc  „m-l  ,  1  \  1  „  \  a  -tine  „m-l 

nxE,  R  dB  +  ~  /  — :V\  n  ■  H  R  dB. 

i-m-1  47r  -t— r*m'.  i-m 

m=°  jb 


Equations  (4. 4a,  b)  and  (4. 5a,  b)  hold  for  all  i  =  0, 1,  2,  . . . ,  however,  the  terms 
i-1  i 

t  and  /*  are  identically  zero  when  i  =  0.  Observe  that  and  are  ex- 
m=0  m=l 


pressed  in  terms  of  the  incident  field  and  preceding  terms  in  the  series  for  E  and 

H,  i.e. ,  E  ,  E, ,  . . .  E .  , ,  H  ,  H. ,  . .  „H  , .  Thus  if  we  consider  the  problem  of 
’  o  l  i-1  0  1  i-1 

finding  E^  and  H^,  assuming  that  the  preceding  terms  have  already  been  determined, 

then  F„  and  G„  are  known  functions, 
i  ( 

The  approach,  ours  as  well  as  Stevenson's,  is  to  show  that  the  unknown 
terms  in  (4.4a)  and  (4.4b)  are  gradients  of  exterior  potential  functions  which  may 
be  determined  as  solutions  of  standard  potential  problems.  Stevenson  went  to  con¬ 
siderable  effort  and  complication  to  formulate  these  problems.  The  method  of  this 
section,  though  still  complicated,  hopefully  represents  a  simplification.  In  any 
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event,  the  present  procedure  for  finding  and  H^,  or  at  least  defining  them  in 
terms  of  solutions  of  potential  problems  is  based  on  the  integral  relation  (4.4a)  and 
(4.4b)  and  does  not  require,  as  Stevenson  does,  first  finding  particular  solutions  of 
Maxwell’s  equations. 

Consider  first  the  task  of  determining  E^(  r).  We  observe,  and  this  is  the 
essence  of  the  approach,  that  the  unknown  term  on  the  right  hand  side  of  equation 
(4.4a)  is  itself  the  gradient  of  an  exterior  potential  function,  a  single  layer  distri¬ 
bution  of  density  n  •  E^.  It  is  possible  to  formulate  a  boundary  value  problem  for 
this  term  as  follows.  Let  0^  denote  the  unknown  potential,  i.e. 


(4.6) 


Then 


E,  “  F,+ 


(4.7) 


where  F„  is  known  and 

i 


V  =  0 


r  in  V 


0^  regular  at  infinity  in  the  sense  of  Kellog 


nx  V0 


a  ,^inc  ri  , 
=  -  nx  (Ef  +FJ) 


r=r 


B 


r-r 


B 


This  is  rot  quite  a  standard  Dirichlet  potential  problem  in  that  the  boundary  condi¬ 
tion  as  given  may  be  shown  to  specify  the  function  0^  on  the  boundary  to  within  a 
constant.  That  is,  specifying  nxV0  on  B  is  equivalent  to  specifying  0  +c  on  B 

l  X 

where  c  is  constant  but  unknown.  This  constant  is  evaluated  by  solving  the  poten¬ 
tial  problem  with  the  ambiguous  boundary  condition,  constructing  the  corresponding 
E^  with  equation  (4.7)  and  then  imposing  the  requirement 
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n- E  HB  =  0 
! 


B 


The  procedure  for  finding  is  thus  seen  to  be  reasonably  straightforward 
once  we  observe  that  the  unknown  part  of  is  the  gradient  of  an  exterior  potential. 
This  observation  spares  us  much  of  the  complication  of  Stevenson’s  approach. 

The  determination  of  requires  more  work  since  it  is  not  obvious  that  the 
unknown  term  on  the  right  hand  side  of  (4.4b)  is  the  gradient  of  an  exterior  potential 
except  when  i  =0.  In  fact  it  may  be  shown  that  when  l  f  0,  this  term  is  definitely 
not  the  gradient  of  an  exterior  potential  function.  Nevertheless  it  is  possible  to 
retain  some  of  the  simplicity  inherent  in  the  determination  of  E^  by  adding  a  known 
function  to  the  unknown  term  such  that  the  sum  is  the  gradient  of  an  exterior  poten¬ 
tial  function.  The  determination  of  the  function  we  must  add  again  requires  the 
solution  of  a  potential  problem . 


Thus  we  introduces  function  g^  as  yet  unspecified,  into  equation  (4.4b),  ob¬ 


taining 


H/r)  *Gi(r)-g^(r)+^  Vx 


nxH, 


R 


dB  +  g|(r)  . 


(4.9) 


'B 


It  is  well  known  and/or  easily  verified  that  a  condition  sufficient  to  guarantee  that 
a  vector  be  the  gradient  of  a  scalar  is  that  the  curl  of  the  vector  vanish,  i.e. 


Vx  A  =  0 


A  =  Vtf/ 


(4.10) 


Thus  a  condition  sufficient  to  guarantee  that  we  may  writ 

1 


nxH, 


47T 


Vx  \  -j-  dB  +  g^r)  =  V^ 
'B 


(4.11) 


is 
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VxVx 


A  g 

nxH. 


dB  +  VXgi  =  0 


(4.12) 


Since  g^  is  as  yet  unspecified,  we  use  (4.12)  as  an  equation  for  and  seek  a  par¬ 
ticular  solution  in  terms  of  known  functions,  i.e.,  terms  not  involving  H  or  E  . 

l  l 

Since  appears,  we  must  first  put  (4. 12)  in  suitable  form.  First  we  use  the  vec¬ 
tor  identity  Vx  Vx  *  V(V  *  -V2  together  with  the  fact  that  for  r  in  V, 


V2  -vr1""-  =  0 

R(r,  rfi) 


to  rewrite  (4.12)  as 


Vxg#(r)  =  VV 

i  47T 


nxH. 


Recall  that  R  is  symmetric  in  r  and  r_  (eqn.  2.1)  and  V—  =  - 

B  R 

operates  on  r  and  Vg  on  r^,  hence  (4.13)  may  be  written 


(4.13) 


VB  ~  ,  where  V 


Vxg/r)  =  VRi  nxH  dB, 


4t r  B  R 

^B 


(4.14) 


or,  on  employing  the  properties  of  the  scalar  triple  product. 


V  x  g^(  r )  = 


^  v \  a-yBRxVB 


(4.15) 


This  we  rewrite  as 


Vxg^(r) 


1  _  \  A  f 

~  V  \  n  •  (V_x— — 
47r  \  B  R 


5  VHi)dB 


(4.16) 
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and,  since  Stokes'  theorem  implies 


(4.17) 


we  have 


Vxgl(r)  =  iv  \  |-Vx§dB 


(4.18) 


But  Maxwell's  equations  (3.10)  imply  that 


VW  =  -Ei-1(?B)  ,>0 


(4.19) 


i  =  1 


hence  we  have,  finally. 


C  n  •  E 

Vxg/?)  -  V  \  dB  . 


f  >0 


(4.20) 


1-0 


We  have  thus  succeeded  in  rewriting  (4.12)  in  terms  of  known  functions  since  we 
have  assumed  that  ^  is  known.  Now  we  want  a  particular  solution  of  (4.  20). 
Clearly  when  i  =  0,  g  =  0  is  a  solution.  When  i  >  0,  we  employ  Stevenson's 
method  for  producing  particular  solutions  of  the  equation  Vx  F  =  f.  First  of  all 
note  that  the  right  hand  side  of  (4.20)  is  the  gradient  of  an  exterior  potential  function 
(single  layer  distribution).  Thus,  introducing  the  notation 
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(4.20)  may  be  written 


Vxg^r) 


reV 


(4.22) 


Stevenson  has  shown  that  necessary  and  sufficient  conditions  for  (4.22)  to  have  a 
solution  are  j(3. 18),  (3.19)] 


V-  VuJ  =  0 


re  V 


(4.23) 


dB  =  0 


(4.24) 


0 

The  first  condition,  (4.23),  is  clearly  satisfied  since,  as  noted,  u^  is  a  potential 
function.  To  show  that  (4.24)  is  also  satisfied  we  use  Gauss’  theorem  to  write 


\  SdB  =  \  s'Vui 

Jr  Jr 


dB  =  - 


\  Vu®dv  + 
1  ‘ 


a  e 
r  •  Vu^  dB 


(4.25) 


where  B^  denotes  a  large  sphere  whose  radius  approaches  infinity.  The  volume 
integral  term  vanishes  by  virtue  of  (4.23)  and  the  surface  integral  over  B  will 

OD 

also  vanish  if 


_  e  ,  2. 

Vu^  =  o(l/r  ) 


as  r— ►go 


(4.26) 


That  (4.26)  is  satisfied  may  be  seen  by  examining  the  structure  of  uj  exhibited  in 
the  defining  equation  (4.21).  Thus 


VuJ(r) 


-  —  V 

4vr 


R(r,  r0) 


(4.27) 


or,  for  r  >  max  r  , 
x3 
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m 

Pm<COS7)dB  <4-28) 

r 

where  cos y  =  cos0cos0_+sin0sin6  cos(0-fL)  and  P  is  a  Legendre  polynomial 

B  15  15  HI 

of  order  m. 

The  m  =  0  term  vanishes  by  virtue  of  (4.3), 


n-  E^dB  =  0 


'-B 


hence 


i  =  0,1.2, . 


co 


Vu  eAi)  = 
l 


4tt 


E 


m=l  r 


m+1 


n-E  (r  )r™P  (cosy)dB 
i- 1  n  u  m 


(4.29) 


'B 


From  (4.29)  it  is  clear  that 


Vu^(r)  =  OO/r^)  or  o(l/r2)  as  r— >oo 


(4.30) 


hence  (4.26)  holds  which  in  turn  means  that  (4.24)  is  valid.  Thus  we  have  estab¬ 
lished  that  equation  (4.22)  has  a  solution.  Furthermore  we  have  shown,  in  the 

3 

process,  that  the  right  hand  side  of  (4.22)  is  0(l/r  )  at  infinity.  Now  we  use  Stev¬ 
enson's  solution  to  this  problem. 

We  define  an  interior  potential  function,  u*(r),  r  interior  to  B,  as  follows. 


V2u*(r  )  -  0 


r  interior  to  B 


n-  VuJ(r) 


r=r 


-  n  •  Vu^(  r ) 


B 


(4.31) 


(4.32) 


r=r 


B 


This  is  a  standard  interior  Neumann  problem  for  u^(  r )  and  has  a  solution  provided 
that 
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1  lei 

Now  using  the  facts  that  V—  =  -V  and  u.  and  u.  are  both  potential  functions  we 

R  v  R  f  i 

obtain 

.  \  V  mB(t  )  C  V  u.(r  ) 

Vxg{  r )  =  ^u(  r  >  -  —  V  V  •  - - -  dv  -  —  V  \  V  •  - - —  dv. 

f  Jf  An  ]  v  R  An  j  v  R 

'v  Jy. 

1  (4.38) 

Now  we  use  Gauss'  theorem,  taking  care  of  the  signs  of  the  normals  (n  on  B  is  al¬ 
ways  directed  into  V,  the  exterior)  to  obtain 

Vx|l(  ?>  -  V»*m  ♦  £  v  \  8  .  [vbu^b>-  VBu‘(rB)]dB  . 


(4.39) 


Actually  there  is  another  surface  integral  term  over  a  large  sphere  at  infinity  but 
this  vanishes  by  virtue  of  (4.30).  The  integral  in  (4.39)  vanishes  because  of  the 
boundary  condition  (4.32)  thus  verifying  that 

Vxg^(r)  =  Vu®(r)  ,  reV  . 

We  may  cast  g  (r)  in  slightly  more  convenient  form  as  follows.  Again 
1  1 

using  the  fact  that  V—  =  -V  —  (4.34)  becomes 

K  V  K 


„  1  _  e.A  ,  ,  1 

V  —  x  V  u  (r  )  dv  -  — — 
v  R  v  I  v  An 


V  ;xV  u*(r  ) dv 
vR  v  I  v 


(4.40) 


or  since  curl  of  the  gradient  is  identically  zero, 


47T 


V  x 
v 


V  u^(r  ) 
v/v 

R 


dv  - 


An 


V  x 
v 


V  u|(r  ) 

V  i  V 

R 


dv  . 


V. 

l 


(4.41) 
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Now  employing  a  famous,  but  apparently  nameless  theorem  of  vector  analy¬ 
sis. 


VxAdv  =  \  nxAds  (S  encloses  V  and  n  is  out  of  S), 


V 

equation  (4.41)  becomes 


l(<?> 


1 

47 r 


nx 


i  B  B 


R 


dB 


f  >  0  (4.42) 


B 


An  alternate  form  of  (4.42)  is  found  to  be 


g/r)  = 


4jt 


Vx 


R 


dB 


(4.43) 


B 


6 

Again  the  behavior  of  Vu^  at  infinity,  (4.30),  causes  a  similar  integral  over  a 
large  sphere  to  vanish. 

0  i 

In  this  form  it  is  clear  that  the  tangential  components  of  Vuf  and  Vuf  on  B 

must  be  unequal  if  is  to  be  different  from  zero.  In  fact  they  are  necessarily  dis- 

0  i 

continuous.  Since  u^  is  a  potential  function  regular  exterior  to  B,  uf  is  a  poten¬ 
tial  function  regular  interior  to  B  and  their  normal  derivatives  were  defined  to  be 

continuous  at  B,  then  the  tangential  derivatives  cannot  also  be  continuous.  If  so, 

i  6 

u^  would  be  a  continuation  into  the  interior  of  B  of  u  .  The  resulting  function 

would  be  a  potential  function  regular  everywhere  in  space  and  therefore  v/ould  nec- 

0 

essarily  be  zero.  But  (see  eqn.  4.21)  is  not  identically  zero. 

We  have  thus  determined  a  particular  g^  such  that  (4.12)  is  satisfied.  This 
in  turn  guarantees  that  equation  (4. 11)  holds,  that  is,  with  the  g  we  have  found  we 
may  write 
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dB  +  =  V\pt  . 


With  equation  (4.34)  or  (4.43)  which  expresses  as  a  curl,  it  follows  upon  taking 
the  divergence  of  (4.  il)  that 

V2^  =  0 . 


On  expanding  l/R  in  (4.43)  it  follows  that 


g^r)  =  0(l/r  )  as 


r— >oo 


A  * 

nxH 


Also 


Vx 


R 


i  2 

-  dB  =  0(l/r  )  as  r— ♦go 


2 

hence  0^  is  regular  in  the  sense  of  Kellog.  (Actually  V0^  =  0(l/r  )  does  not  imply 
completely  that  0^  is  regular.  There  may  be  an  additive  constant  which  would  im¬ 
ply  jr0^|  is  not  bounded.  Since  we  are  interested  in  V0^.  which  removes  this 
constant  anyway,  we  may  choose  it  as  zero  to  begin  with  and  take  0^  to  be  regular.) 

With  equation  (4.9)  and  the  boundary  conditions  (4.2)  we  may  formulate  a 
standard  exterior  Neumann  potential  problem  for  0^,  namely 

V20i  =  0  re  V 

0^  regular  at  infinity  in  the  sense  of  Kellog  (4.44) 


A 

n 


V0! 


«4  ■* 

r=r 


30i 

9n 


B 


r=? 


a  £ine  a.* 

-n  •  H  -  n  •  G  + 
I  £ 


B 


“xVui  - 


R 


i’l 


j 


r=r 


B 


The  solution  of  this  problem  then  determines  H^(  r). 
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To  summarize  the  procedure  we  have  established: 
If  an  electromagnetic  field, 

CD  CD 


Eta°< 


f)  EfC(r)(ik)'  .  Hinc(f)  =  ^Hfc(?)(ik)< 

1=0  1=0 


(4.45) 


is  incident  on  a  smooth  finite  perfectly  conducting  surface  B  in  three  space  then  the 
coefficients  in  the  low  frequency  expansion  of  the  scattered  field 


oo 

E^rXik) 

are  given  by 

E/r)  =  F/{r)  +  V0i  (4.47) 

H^r )  =  Gp)  -  gp)+  (4.48) 


H(r)  = 


/=0 


H^(  r  )(ik) 


(4.46) 


where 
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(we  must  use  \ft  •  E^dB  =  0  to  determine  an  arbitrary  constant  arising 
from  this  form  oT  the  boundary  condition), 

and  0  is  an  exterior  Neumann  potential: 


A  - 0 

^  reg  at  a> 


r  exterior  to  B 


A  „ ,  a  *mc  a  -±  a  -» 

n  •  V0^  =  -n  •  -  n  •  G^+n  •  on  B 


(4.55) 


J 


n  is  the  unit  normal  on  B  always  pointing  into  V,  the  exterior  of  B,  and  R  is  the 
distance  Jr  -  ?BJ  from  a  point  rfi  on  the  surface  (the  integration  variables)  to  a 
field  point  r. 

We  complete  this  section  with  a  brief  discussion  of  the  low  frequency  expan¬ 
sion  of  the  far  field.  Here  we  proceed  exactly  as  in  the  scalar  case.  We  incorpor¬ 
ate  the  facts  that,  for  large  r. 


ikR 

e  e 

'  -i  .I 


R 


ikR 

p  ^  ..  a  e 

V — — —  ~  liar  - - 

R  r 


in  the  Stratton-Chu  integral  representations  of  the  scattered  field,  equations  (3.4a) 
and  (3.4b),  also  employing  the  boundary  conditions  on  the  surface  (3. 12),  obtaining 


E(  r )  ^ 


H(r) 


ikr 

_e _ lk 

4?rr 


ikr .. 
e _ lk 

4rrr 


-ik 


r  •  r 


B 


*  a  ^inc,  ,  a  a  a  * 
-rx(nxE  )  +  nxH-rn-E 


dB 


(4.56) 


B 


-ikr  •  r 


B  Ta  ,a  a  .  ^inc.  a  a  r*incl 
jrx(nxH)  +  nxE  +rn-H  J 


dB  (4.57) 


IB 


Now  we  expand  the  field  quantities,  equations  (4.45)  and  (4.46)  and  the  factor 
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-lkr  •  r. 


and  reorder  the  terms  to  find 


E  r  4ffr  Y7T*  ^ 

1=0  m=0 


(r  •  rfi)  |-rx(nx  E^ 


+  nxH 


-rn  •  E  ~|c 
i  mj 


(4.58) 


«<?>  IJtjS  \  <*V‘m 

<=0  m=0  Jg 

+  nxElnC  +  rn-  H  ““HdB  (4.59) 
m  m  J 

The  1  =  0  term  in  (4. 58)  and  (4.  59)  always  vanishes  the  integrals  \  n  *  E  dB 
r  L  m 

and  \  n  -  H  are  zero  for  all  m,  (equation  4.3)  and,  since  E  and  H  may 
m  o  o 


both  be  written  as  gradients  of  potential  functions,  we  may  use  a  well  known  re¬ 
sult  of  vector  analysis  which  implies  that  ^fix  VjJdB  =  0,  for  B  closed  and  any 

0  to  see  that  \  fixH  dB  and  \  nxEmcdB  also  vanish.  .  Therefore  we  may 

4  0  4  0  J 


rewrite  (4. 58)  and  (4. 59)  as 


E(r) 


k2 

47rr 


'y' ok)' f1, [ <*  - 

u  «+i-*  Vr-r> 


J+l-m 


•  f-^x^xE*110  )  +  ?lxH  -rn*E  ~]dB 
L  m  m  u\J 


(4.60) 
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H(r) 


''V  - 


e^k2 

4  tT 


■f5J<lk)  (i+l-m): 
1=0  m=0 


(r  -  rB) 


f+l-m 


'B 


A  i*  «  \  ,  *  ^  inc  a  ^  -me 
rx(nxH  )  +  nxE  +  rn*H 
m  m 


'JdB  . 

(4.61) 


This  illustrates  a  famous  result  of  Rayleigh:  the  leading  term  in  a  low  frequency 

2 

expansion  of  the  far  field  is  proportional  to  k  .  Stevenson  criticized  this  form  as 
being  inefficient  since  one  apparently  needs  to  determine  f+1  non-vanishing  near 
field  terms  in  order  to  obtain  i  non- vanishing  far  field  terms.  Actually  this  is  not 
completely  true,  as  a  close  examination  of  the  "extra”  near  field  terms  reveals. 
These  are  the  m  =  f+1  terms  in  (4.60)  and  (4.61),  namely 


Ea  .fs  inc.  .  a  —  A  A  -  *"] 

rx(nxEw)  +  nxHm-rn-EwJdB 


(4.62a) 


and 


[?x(axHm)  +  &xE™  +  n-H“°]d 


(4.62b) 


IB 


which  we  rewrite  as 


A 

-rx 


a  ->inc  \  a  a  \  a 

nxE^dB+\  nxH|+1dB-r\  n  •  x  dB 


(4.63a) 


Jr 


IB 


and 


nxK,,  dB  + 
1+1 


a  ^rnc  ,  a 
nxEj,.  dB  +  r 

rr  1 


a  *inc 
n  •  dB 


(4.63b) 


'B 


JB 


IB 


The  terms  involving  the  incident  field  are  effectively  known  since  the  inci- 

ffl 


dent  field  is  given.  Also,  from  equation  (3.14),  \  n  •  dB  =  0  and  the  only 


:ia  are  eueci 

1.14),  \  n 
JB 
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unknown  part  of  these  "extra"  terms  involves 


nxH^,  dB  . 
i+1 


(4.64) 


'B 


With  (4.48),  however,  it  follows  that 


fixH  dB 

i+1 


B 


r\ 

A 

J 


(nxGm-  nxV^ldB 


(4.65) 


But 


l 


nxV^  dB  =  0,  using  a  vector  identify  we  have  employed  before,  and 

rrl 


G  is  given,  (4.50),  in  terms  of  the  first  l  near  field  terms.  The  only  unknown 

jrrl 

part  of  this  "extra"  term  involves  g  which  does  require  the  solution  of  an  inter- 

rrl 

ior  Neumann  problem  see  equations  (4.  51)-(4.53)  .  This  is  considerably  less  than 
requiring  complete  determination  of  E.,  and  H  ,  but  is  still  unsatisfactory. 
Repeated  attempts  to  determine  this  "unknown"  part  without  solving  for  g^  have 
so  far  been  fruitless.  The  alternatives  are  also  less  than  overwhelmingly  desirable] 
Stevenson  provides  a  generalization  of  Rayleigh's  continuat  on  method  whereby  the 
near  field  terms  for  large  r  are  matched  with  multipoles  for  small  k  (thus  defining 
the  multipole  moments)  then  using  the  far  fields  of  the  multipoles.  This  of  course 
involves  expanding  the  near  field  terms  in  spherical  harmonics  which  may  involve 
as  much  labor  as  solving  the  required  interior  Neumann  problem.  Still,  in  princi¬ 
ple,  Stevenson's  method  of  continuing  into  the  far  field  is  preferable  since  it  does 
not  require  the  solution  of  another,  albeit  simple,  problem  in  order  to  obtain  the 
same  number  of  terms  in  a  low  frequency  expansion  of  the  far  field  as  are  available 
in  a  low  frequency  expansion  of  the  near  field.  The  price  is  apparently  requiring 
both  to  be  represented  as  expansions  in  spherical  harmonics 


41 


THE  UNIVERSITY  OF  MICHIGAN 

7133-4-T 


AN  EXAMPLE— SCATTERING  BY  A  SPHERE 

To  illustrate  the  procedure  derived  in  the  previous  section,  we  consider  the 

problem  of  scattering  of  a  linearly  polarized  plane  wave  by  a  sphere.  The  incident 

-*inc 

field  is  taken  to  propagate  down  the  z-axis,  with  E  along  the  x-direetion  (see 
Fig.  2),  i.e. 


-me  a  -lkz 
E  =  1  e 
x 


uj 

9*“  =  A  e'ik2  = 


-inc_  (-z)  * 
E<  -  —  ‘x 


-inc  _  (-z)  a 

Hi - FT  V 


(5.1) 


H  ™ 


EmC 


Figure  2 


/B=a 


We  shall  proceed  to  calculate  the  first  two  terms  in  the  series  for  the  scattered 
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op 

E  ~^>_,<ik/E,  . 


gp 


(5.2) 


by  straightforward  application  of  equations  (4.46)  -  (4. 55). 


5. 1  Zeroth  Order  Terms 


From  (5. 1)  we  see  that 


EtaC  =  t  . 

O  X 


(5.3) 


The  scattered  electric  field  to  this  order  is  (4.47) 


E  =  F  +  V0 
o  o  o 


(5.4) 


and  (Eq.  4.49) 


1  \ 

F  =  -7- Vx  \  — ~  dB 
o  4tt  \  R 


(5.5) 


Equation  (5.5),  written  in  its  entirety  is, 


F  (r)  = 
o 


p27 r  p?r 

”  V  x  \  d0fi  \  dSBa28inSB 


A  A 

rBX‘x 


2  2  r  "i 

r  +a  -  2ar  |cos6cos0B+sin6sin6Beos(0- 0g)| 


(5.6) 


where  the  unit  normal  is  r  =1  sin0  cos0  +i  sin0  sin0  +i  cos0  . 

J3x  is  By  is  is  Z  is 

The  integration  is  carried  out  using  the  well  known  expansion  of  l/R  in 

spherical  harmonics  and  we  find  (using  a  mixture  of  rectangular  and  spherical  unit 

vectors) 
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From  (4.50) 


H  =  G  -rV0  . 
o  o  o 


\A  £inc 

n-  H 

-r~dB 


(5. 14) 


(5.15) 


from  which 


=  4 

o  y 


G  = 

o  3  \  3  / 


(5. 16) 


(5.17) 


ii  (Eq.  4.55)  is  an  exterior  potential  function  with  boundary  values 

CH>J  .  .  a  /  3  \ 

o  a  ,-trnc  a  £  a  a  ,  1  9  /  a  \ 

—  =  -r  •  H  -  r  •  G  =  r  •  l  +  -  —  ~y 

3r  o  o  y  3  3r  \  3  / 


=  —  sin#  sin  0. 

Such  a  function  is  easily  seen  to  be 


(5. 18) 


1  3  ** 

(p  =  -  -  —  sin#sin0 
o  b  2 
r 


which,  with  (5.14)  and  (5.17)  leads  to 


i  /  3  \ 

H  =  -  -  V  (  j  = 
o  2  \  3  J 


a3!  „  3 

— *  +  - 
,  3  2  4  r 

2r  r 


(5.19) 


(5.20) 


The  zeroth  order  results  may  be  rewritten  entirely  in  terms  of  spherical  unit  vec¬ 
tors  as 
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*  tea2  Aa3  a  a3 

E  =  — r-  sin0cos0- 6  —z  cos0cos0+0  —  staff 

0  o  o  0 

r  r  r 


(5. 21 


H  =  •~-8ln£)8in0- 6  cos8sta0-0  cos0 


(5.25 


5.2  First  Order  Terms 

The  next  terms  are  found  using  these  results,  again  following  the  procedure 
of  the  preceding  section.  It  is  to  be  noted  that  even  at  this  stage,  the  calculations 
become  tedious.  With  (4.49)  we  see  that 


V 


J  C  nxSj110  C  f  nxH 

-^Vx\  — — °6B+^  \  ~irdB 


— -  V  \  n  •  E  dB 
4x  \  o 


(5.23 


J  Jk  Jklnp 

All  terms  are  well  defined,  E  and  H  in  (5.21)  and  (5.22)  above  and  E  and 
^inc  o  o  o 

in  (5.1),  namely 


EtaC  =t 

0  X 


>mc  a 
E,  =  -zi 
1  x 


(5.24 


Carrying  out  the  indicated  integrations  we  find 


n  x  z  a  a  a  ,  a  /  ^  ,  ft  \ 

F  = - T  r- ~~  rxi  +  — -  (zi  +Xi  ) 

1  _6  „  2  y  5  x  z 


(5.25 


Now  we  use  (4.54)  to  find  0^.  The  boundary  condition 


a  a  -inc  =*  . 

nxV0  =  -nx(E,  +F  )  , 

1  x  1 


is  seen  to  imply 
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3a2  a2  1 

P!  =  IF  cos0sin2@+c  =  ~Y5  P Acqb6)cqb$  +c 
r=a 


(5.26) 


Here  we  employ  the  definition  of  the  associated  Legendre  functions  given  by  Magnus 
and  Oberhettinger  (1949).  The  exterior  potential  function  taking  on  these  boundary 
values  is  found  to  be 


a5  > 


0  _  -  p  (cost?)  COS  0  +  — 

1  10r3  2  r 

Forming  with  (4.47),  and  applying  the  auxilliary  condition  \ 
which  implies  c  =  0,  we  find 


(5.2?) 


dB  =  0 


3  ,  /5 

1  =  +iv(i4 

1  2r2  y  2  V,5 


(5.28) 


Proceeding  to  the  determination  of  H  we  see  that  (Eq.  4.48) 


*imQr*i+v*i- 


(5.29) 


With  (4.50) 


*•  1  \  A  1 

G  =  ~  Vx  \  r  xH  dB  + — 
1  47T  \  B  o  4tt 


r  xEmC  ,  *  . 

- -  dB  +  ~  V  — J.  ^ 

B  47 r  R 


1  a  -*inc 

+  ¥V  VH,  «■  (5-30) 


All  the  quantities  involved  have  already  been  defined  in  (5. 1)  and  (5.22).  Carrying 
out  the  integration  yields 


n  a  *  a 

o  =  — :«> 

3r“  ’ 


‘x  +  v(fs  y*)  • 

\5r  / 


(5.31) 
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6  i 

To  determine  g^  we  must  first  find  and  u^.  From  (4.52) 


1  f 

e  1  \  o 

U1  ”  4tt  j  R 


dE 


(5.32) 


B 


which,  with  (5,21),  may  be  evaluated  as 


U1  ~  f  ^2  P^CO8^CO80 

r 


(5.33) 


Following  (4. 53)  we  determine  the  interior  potential  function  whose  normal  deriva- 

0 

tive  matches  that  of  on  the  boundary.  Here  we  assume  a  series  of  the  form 

a  rDP  (cosy)  and  determine  the  a  using  the  boundary  condition.  In  the  pre- 
n  n  d 


sent  case,  this  is  easily  seen  to  be 


i  4  1 

u  =  - -r  P  (cos0)cos0  . 

X  0  1 


(5.34) 


With  (4. 51)  we  see  that 


8,  = 


nxVfi(u®  -u* ) 
~~  R 


dB 


(5.35) 


and  using  the  expressions  (5.33)  and  (5.34'  this  may  be  found  explicitly  as 

3 

•v  2  a  a  j 
g,  =  -  “  —  rxi  . 


'1 


3  2 

r 


(5.36)1 


Now  we  proceed  using  (4. 55)  to  find  an  exterior  potential  function,  0  ,  satisfying 
the  boundary  condition 


* \ 

Bn 


■A-h*0  A 


A 


r=a 


p -n.G1+n-g1 


(5. 37)1 
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With  (5. 1),  (5.31)  and  (5. 36),  this  boundary  condition  becomes 


"V'l  2  2a  2 

T  =  -  -  asin0cos0sin0  =  —  PJcos0)  sin0. 
dr  o  i j  Z 

r=a 


(5.38) 


and  the  solution  is  found  to  be 


p2(cos0)sin0  =  — ^ 
45r  15 r 


(5.39) 


Substituting  (5.31),  (5.35)  and  (5.39)  in  (5.29)  we  obtain 


3 

a  a  * 
~t:  rxi 

2  i 

r 


t  +iv(^ 

x  3  V  5 


(5.40) 


The  first  order  terms  may  be  written  entirely  in  spherical  coordinates  as 


5  /  5  3 

E  =  -  —  ~  sin0cos0cos0r  +  f  cos  20+ 

1  r  '2r  2jl' 


/  3  5 

-  cos0sin0  ( 

\2r  2r 


(5.41) 


5  /  5  3\ 

=  -  ~  sin0cos0sin0  r  +  (-—  cos  20+  — ^  )  sin0^ 
r  '3r  r  / 


+  COS0COS01  —  +  — - 
\r  3r 


(5.42) 


These  results  for  the  first  two  terms  in  the  low  frequency  expansion  may  be  shown 
to  be  in  complete  agreement  with  comparable  expressions  derived  from  the  standarc 


Mie  series. 
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